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Thick gas discs in faint dwarf galaxies 
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ABSTRACT 

We determine the intrinsic axial ratio distribution of the gas discs of extremely faint 
Mb > —14.5 dwarf irregular galaxies. We start with the measured (beam corrected) distri- 
bution of apparent axial ratios in the HI 21cm images of dwarf irregular galaxies observed as 
part of the Faint Irregular Galaxy GMRT Survey (FIGGS). Assuming that the discs can be 
approximated as oblate spheroids, the intrinsic axial ratio distribution can be obtained from 
the observed apparent axial ratio distribution. We use a couple of methods to do this, and our 
final results are based on using Lucy's deconvolution algorithm. This method is constrained to 
produce physically plausible distributions, and also has the added advantage of allowing for 
observational errors to be accounted for. While one might a priori expect that gas discs would 
be thin (because collisions between gas clouds would cause them to quickly settle down to a 
thin disc), we find that the HI discs of faint dwarf irregulars are quite thick, with mean axial 
ratio < q >^ 0.6. While this is substantially larger than the typical value of ~ 0.2 for the 
stellar discs of large spiral galaxies, it is consistent with the much larger ratio of velocity 
dispersion to rotational velocity {a/v c ) in dwarf galaxy HI discs as compared to that in spiral 
galaxies. Our findings have implications for studies of the mass distribution and the Tully - 
Fisher relation for faint dwarf irregular galaxies, where it is often assumed that the gas is in a 
thin disc. 
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1 INTRODUCTION 

The intrinsic shape of galaxies is interesting from a variety of per- 
spectives. For a given galaxy the shape should be consistent with 
the dynamical model of the galaxy, while, for a sample of galax- 
ies, one would expect that a correct evolutionary model should 
be able to reproduce the observed distribution of shapes. It is 
generally assumed th at disc galaxi e s can be approximated to be 
oblate s pheroid s (e.g. [Hubble! l ll926t) : ISandage. Freeman. & Stokes] 
(1970); Rvden (2006)). If one further assumes that the galaxies 
have a well defined mean axial ratio (qo), then the observed ax- 
ial ratio can be used to determine the inclination of the disc. In 
turn, the inclination is a crucial input in dynamical modeling (e.g. 
for the mass and structure of the dark matter halo), studying the 
Tully-Fisher relation etc. 

The observed shape of a galaxy differs from the intrin- 
sic shape because of projection effects. If one has a sample of 
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galaxies drawn from a population with a well defined intrin- 
sic axial ratio distribution and with random orientations with 
respect to the earth, then one can determine the distribution 
of intrinsic a xial ratios from the observed axial ratio distribu- 
tion (for e.g. lNoerdlingerlll979l: iBinnev & de Vaucouleurslll98ll: 
lLambas, Maddox & Lovedavlll992l : lRvdenll2006h . It is worth not- 
ing that most of these studies have focused on large galaxies, 
and that there have been relatively few that focused on dwarfs. 
For bright spiral galaxies, qo i s often taken to be ~ .2 (see 
e.g. lHavnes & Giovanelli|[l9p ; IVerheiien & SancisfcOOlh . It has 
also long been appreciat ed that the axial ratio is a function of Hub- 
ble ty pe. For example, iHeidmann. Heidmann. & de Vaucouleursl 
( 1972) found that while discs get thinner as one goes from galax- 
ies of morphological type Sa to Sd, there is a rapid increase in disc 
thickn ess as one goes from Sd to dwarf Irregu lar galaxies. Sim- 
ilarly, jstavelev-Smith, Davies & Kinmar] i 19921) found that dwarf 
galaxies from the UGC catalog have qo ~ 0.5. Axial ratio is also 
a func tion of the wavelength of observation. For example, Ryden 
(2006) showed that older populations as traced by redder stars have 
thicker ratios than the corresponding B band disc. But all of these 
studies refer to the stellar discs of the galaxies. Due to collisions 
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Figure 1. Histograms of the measured apparent axial ratios (after correc- 
tion for the finite beam size). There are a total of 15 bins, the error bars 
were determined assuming Poisson statistics. The shaded area represents 
the 80% confidence limit histograms generated by bootstrap re-sampling of 
the data (see text for details). The dashed line shows the best fit (i.e. using 
the weighted least square) polynomial to the binned data. The solid black 
line is the reconstructed distribution (<£(p)) given by Lucy deconvolution 
(see text for details). All values have been normalized so that the area under 
the curves and histogram is unity. 



between gas clouds, one would expect that the gas discs would be 
intrinsically quite thin, for example, for our Galaxy, the scale height 
of the middle disc is only ~ 300 pc. In extremely gas rich dwarf 
galaxies one might then expect that the gas discs are relatively thin, 
even though the stellar disc may have a large axial ratio. 

We present here a study of the axial ratio distribution of 
the HI discs of extremely faint dwarf irregular galaxies. We se- 
lect a sample of galaxies for which HI synthesis observations are 
available from the F aint Irregular Galaxy GMRT Survey (FIGGS, 
iBegum et all (2008)). The sample selection and analysis methods 
used are presented in section [2] while the results are discussed in 
section [3] To the best of our knowledge, this is the first such study 
of the intrinsic shape of the HI discs of galaxies. 



2 SAMPLE SELECTION, ANALYSIS AND RESULTS 

The FIGGS sample was selected to satisfy the following criteria: 
(1) absolute blue magnitude Mb > —14.5, (2) integrated HI flux 
> 1 Jy km/s and (3) optical major axis > 1 . There are total of 62 
galaxies in the FIGGS sample, for all of which we have axis ratios 
available. We discuss the implications of the selection criteria on 
our derived axis ratio in Sec. [3] here we merely note that any bias 
in the sample is likely to be towards having a preference for edge 
on systems. The FIGGS survey presents HI images at a range of an- 
gular resolutions. Since we are interested in the overall shape of the 
disc, we measured the axis ratio of the disc using the coarsest reso- 
lution HI maps. This was done because we wanted to avoid the fine 
structure which shows up at higher resolutions, and instead trace 
the complete spatial spread of HI which requires observing the less 
dense gas using coarser resolution. For all of these galaxies, ellip- 
tical isophotes were fit to the HI image using the STSDAS package 
in IRAF. The axis ratio that we use in this paper was measured for 
the isophote with mean column density 10 19 cm -2 for all galax- 



ies. The measured axial ratio was corrected for the finite beam size 
before being used. 

For a sample of randomly oriented axi-symmetric oblate 
spheroids, with intrinsic axial ratio distribution ip(q), where the 
axial ratio q is defined as q — c/a, with c being the short axis 
("thickness") of the spheroid and a the long axis ("diameter"), the 
observed distribution of axis r atios <j>(p), wh ere p is the observed 
axial ratio is given by (see e.g. lBinnevI jl978h ) 



Tp(q)dq 
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What is available to us is 4>(p) the observed axial ratio distri- 
bution; in order to derive the intrinsic axial ratio distribution ip(q) , 
we h ave to invert Eqn.[T] 

iFall & Frenkl dl983f) provide an integration formula (Eqn. (7) 
in their paper) in the case that cj>(p) is a polynomial function. This 
allows for a direct inversion of Eqn. [T] The binned data (Fig. [TJ 
were hence fit by a polynomial, using a weighted least squares fit. 
The errors were assumed to be Poisson distributed (with the er- 
ror for bins with zero counts taken to be that for a single count). 
The errors were also estimated using bootstrap resampling. 10000 
different realizations of the data were constructed using boot strap 
resampling, and the 80% confidence limit histograms determined 
in this way are also shown in Fig. Q] As can be seen, the two er- 
ror estimates are in good agreement. The direct inversion of the 
best fit polynomial is shown in Fig. [2] It has a broad peak around 
qo ~ 0.6, but also falls to negative values beyond q ~ 0.8, which 
is unphysical. 

While there may exist an acceptable polynomial fit to 
<j>(p) which does not produce negatives on inversion, we do not in- 
vestigate this possibility directly. Instead we cho ose to try to invert 
Eqn. Q] using Lucy's deconvolution algorithm dLucvlll974l) . This 
approach has the added advantage of also being able to account for 
measurement errors. The Lucy deconvolution method is an iterative 
algorithm for estimating the intrinsic frequency distribution from 
an observed distribution, subject to the constraints that the deduced 
distribution should be normalized as well as positive for all values 
of the intrinsic quantity (q in our case). One starts with an initial 
guess initial guess for the intrinsic distribution (tp(q)), and uses a 
kernel (k(p\q), defined by Eqn. [2] below) to iteratively find better 
approximations to both the apparent distribution and the intrinsic 
distribution. 

If we rewrite Eqn. [T] as: 



(p{p) = J ip(q)k(p\q)dq 
then the rth loop of the iterative algorithm will be, 

#-(p) = [ iP r (q)k(p\q)dq 
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Eqn.|4]neglects the mea surement errors in determining t he ob- 
served axial ratio. Following Binnev & de Vauco uleursl (1981), wc 
can account for these errors by assuming the following form for the 
probability E(p\p')Sp, that a galaxy of actual apparent axial ratio 
p' is instead measured to have an axial ratio p: 



E(p\p') = 
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Figure 2. The frequency distribution of intrinsic axial ratios. The dashed 
line indicates i/j(q) as obtained from the direct inversion of the polynomial 
fit to the binned data. The solid line is the final ip(q) obtained after Lucy 
deconvolution (see text for details). The shaded area represents the region 
between the ip(q)s obtained by applying Lucy deconvolution to the polyno- 
mial fits to the 80% confidence interval histograms (see text for details). 



In our case, we estimate the measurement errors by fitting a 
straight line to the distribution of errors in the axial ratio p as re- 
ported by the isophote fitting package. This gives: 



cr(p') = 0.029p' + 0.0006 
and projection kernel hence becomes, 



K(p\q) 



E(p\p')k{p'\q)dp 



(6) 



(7) 



The algorithm now involves using K(p\q) instead of k(p\q) in 
equations [3]and[4] for which the range of integration now becomes 
to 1, since now K(p\q) 7^ for p < q. One would expect that 
in successive iterations, the approximation to the observed 
distribution (f)(p) improve, and this can be used to decide when to 
terminate the algorithm. 

Two different initial guesses for ip(q) 0/>°(?)) were tried, the 
first was peaked around q ~ 0.6, the second was constant for all q. 
These two initial guesses thus cover two extremes of possible tp(q). 
For both initial guesses, ip(q) quickly converges to the form shown 
in Fig. [2] To determine the indicative errors in the determination of 
ip{q) , the 80% confidence limit histograms shown in Fig. [TJ were 
fit with polynomials and Lucy deconvolution was applied to these 
polynomials. The resulting intrinsic distribution functions are also 
shown in Fig. [2] and the shaded area hence represents the ~ 80% 
confidence interval for tp(q). The final <&(p) (i.e. the final approxi- 
mation to the observed distribution function) is shown in Fig.[TJ and 
as it can be seen, it fits with the observed distribution within the er- 
ror bars. The ip(q) shown in Fig.|2]has < q > = 0.57, standard 
deviation a = 0.164, and skewness, 71 = q s /a s — —0.62. 



3 DISCUSSION 

Eqn. [T] (and hence results obtained from inverting it also) applies 
only for a sample of randomly oriented oblate spheroids. The basic 



assumption hence is that the galaxy sample we are working with 
has an unbiased distribution of inclination angles. The selection 
criteria for the FIGGS sample (from which the current sample is 
drawn) includes a requirement that the optical major axis of the 
gala xy be > 1 . Dwarf galaxies are genera lly dust poor (for eg. 
see IWalter et al.ll2007l ; lGalametz et al. 2009). and hence to a good 
approximation have optically thin discs. Highly inclined galaxies 
will hence be over represented in a diameter limited sample, i.e. 
our sample is biased towards edge on discs. This means that the 
true mean intrinsic axial ratio would be even larger than what we 
have estimated above. It is worth noting that as the intrinsic axial 
ratio gets closer to 1.0, the magnitude of this bias decreases, and 
hence the bias in our estimate should not be large. 

The mean intrinsic axial ratio that we obtain, viz. < q >~ 
0.57 is substantially larger than the value of 0.2 usually adopted 
for the stellar discs of large spiral galaxies. The value from 
our sample is more than twice as large as older measurements 
of < q > in stellar discs of Magellanic irre gular galaxies 
by Heidmann, Heidm ann. & de V aucouleursj J 1972b (ranging from 
0.20 to 0.24). Consistent with this, our sample contains no very 
flat galaxies. In fact, as can be seen from Fig. [T] there are no 
galaxies with apparent axial ratio < 0.2 in the sample. Further, 
if we look at the distribution of the observed axial ratio of differ- 
ent classes of galax ies in the Automated Photographic Measuring 
(APM) survey as inlLambas, Ma ddox & Lovedavl dl992l) . then our 
histogram resembles those for ellipticals and SOs more closely than 
that for spirals. This is another qualitative indication that the un- 
derlying intrinsic distribution of axial ratios has a higher mean than 
is typical for spirals. I nterestingly, the value of < q > we obtai n 
matches well with whatlstavelev-Smith. Davies & Kinmanl ( fl992r) ; 
Binggeli & Popescul dl995l) derived for the stellar discs in dwarfs. 

The thickness of the gas discs of dwarf galaxies is contrary 
to what one might have naively expected for a gas disc, since, 
in general, collisions between gas clouds should cause them to 
quickly settle into a thin disc. However, this large axial ratio is 
probably consistent with the large gas dispersion in comparison 
to t he rotational velocity observed in dwar f galaxies. For exam- 
ple. lKaufmann. Wheeler and Bullockl (2007) did particle hydrody- 
namic simulations to show that dwarf galaxies with rotational ve- 
locities ~40 km s _1 did not originate as thin discs but thick sys- 
tems. This still leaves o pen the ques t ion of where the large velocity 
dispersion comes from. iDutta et all d2009h find good evidence for 
a scale free power spectrum of HI fluctuations in dwarf galaxies, 
consisent with what would be expected from a turbulent medium. 
Interestingly, they also find that that the dwarfs must have relatively 
thick gas discs, similar to the conclusions reached here. Assuming 
that the origin of the velocity dispersion is turbulent motions in the 
ISM, the timescale for dissip ation is given by r ~ L/v tur t ~ L/a 
dShu. Adams. & Lizano|[l987i e.g .). The total turbulent energy is 
Eturb ~ 1/2Mhio~ 2 - For our sample galaxies, the typical HI mass 
is ~ 3 x 10 7 Mq, while the length scale L ~ lkpc. The rate of tur- 
bulent energy dissipation is hence ~ 10 44 erg/ yr. On the other hand, 
the s tar formation rate is ~ 10~ 3 M0/yr dRovchowdhurv et al.l 
2009), for which t he expected su pernova rate for a Salpeter IMF 
is ~ 7 x 10 _6 /yr dBinnevlfeoOll) . Assuming that each supernova 
explosion deposits ~ 10 51 ergs of energy into the ISM, the energy 
input from star formation is ~ 10 46 ergs/yr, more than sufficient to 
balance the turbulent energy loss. Thus, star formation driven tur- 
bulence in the ISM is a plausible cause for the thick gas discs that 
we observe. 

We have assumed through out that the HI discs of dwarf galax- 
ies can be approximated as oblate spheroids. On the other hand, 
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Figure 3. Reconstructed axial ratio distribution (<I>(p)) assuming the gas 
discs to be prolate spheroids, is shown as boxes, which follows the best 
fit polynomial to the binned data (solid line) closely. Binned data with er- 
rorbars and the reconstructed distribution &(p) obtained using the oblate 
spheroid assumption (dashed line) is shown for comparison. All values have 
been normalized to unity. 



in Sec. [2] we saw that the inversion based on the best fit polyno- 
mial to the observed histogram of axial ratios actually gives un- 
physical results, i.e. that the derived intrinsic axial ratio d i stribu - 
tion ip{q) becomes negative. lLam bas. Maddox & Lovedavl d 19921) 
found a similar pattern for the axial ratio distribution of spiral 
and SO galaxies in the APM catalog, and hence relaxed the as- 
sumption that the discs are oblate spheroids. Adequate fits to their 
data could be ob tained assumin g that the galaxies have a triaxial 
shape. Similarly, Ryden ( 2006) from a study of large galaxies in 
the 2MASS catalog concluded that spiral galaxies are mildly tri- 
axial. The fact that the polynomial approximation gives unphysi- 
cal results for our sample also suggests that triaxial models may 
provide a better fit. On the other hand, the Lucy deconvolution 
gives a physically plausible (as indeed it is designed to) intrinsic 
axial ratio distribution, that fits the observed data within the error 
bars. It is worth noting that distribution found by Lucy deconvo- 
lution is in excellent agreement with that found by direct inver- 
sion of the polynomial fit for the entire range for which the lat- 
ter is > 0. Interestingly, if we assume the gas discs to be pro- 
late instead of oblate spheroids, Lucy deconvolution produces an 
equally acceptable <j>(p) as can be seen from Fig. [3] The ip(q) ob- 
tained from Lucy deconvolution with a prolate spheroid assump- 
tion (see Fig.0 indicates that such gas discs should be more cylin- 
drical than spheric al. However, the observed kinematics of gas in 



these ga l axies (see iBegum et alj|2003l; Begum & Chengaiml 



2004affbl: IBegum Chengalur & Karachentsevl |200^ ~ Begum et 



2003, 
et al.l 



2005 1120061" I2008L etc.) shows that the gas has a significant rota 
tional support. As discussed earlier, a gravitationally bound rotat- 
ing disc of gas forms an oblate and not a prolate spheroid. Thus, 
although from the axial ratio data alone, one cannot distinguish be- 
tween a prolate and oblate shapes, in conjuction with the kinemat- 
ical data, it is clear that the gas is distributed in the form of a thick 
disc. 

In summary, we find that the gas discs of dwarf galaxies are 
relatively thick, in sharp contrast to the gas discs of spiral galaxies. 
This has implications both for the internal dynamics of the gas, as 
well as for studies of the mass distribution Tully - Fisher relation 
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Figure 4. The frequency distribution of intrinsic axial ratios obtained after 
Lucy deconvolution. The dashed line indicates ip(q) obtained assuming 
the gas discs to be prolate spheroids, while the solid line shows the earlier 
ip(q) obtained with the oblate spheroid assumption. 



in faint dwarfs, in which it is generally assumed that the gas is in a 
thin disc. 
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